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Abstract. In this note we consider /c-regular multigraphs, where the possible edge multiplicities are con¬ 
trolled. These structures are considered in a question recently posed by Brendan McKay. We express the 
generating functions using the scalar product of symmetric functions, and consequently give conditions on 
when the classes are D-finite. We appeal to symmetric species results of Mendez to write the expressions in 
a systematic way. 


1. Introduction 

The asymptotic enumeration of regular graphs is a compelling topic that has appeared in many forms in 
combinatorics over the past half century. There are several approaches, and each has its own conditions, 
and results. In this note we revisit the symmetric function approach, introduced by Goulden, Jackson and 
Reilly [7] , generalized by Gessel 0, and automated by Chyzak, Mishna and Salvy [3]. The goal of this work 
is to give insight on a recent related problem posed by McKay: 

Problem (McKay [12]). Let 3, and 3* be subsets of the non-negative integers, and let d = d{n) = (di ,..., d n ) 
be a vector of non-negative integers. Let M(n,3,3*) be the number of symmetric matrices whose diagonal 
entries are drawn from 3* and off-diagonal entries from 3 , whose row sums are d \,..., d n . As usual in graph 
theory, entries on the diagonal are counted twice. We are interested in the asymptotic value of M(n,3,3*) 
in the sparse case, where the row sums do not grow very quickly with n. 

We consider this problem in the case that the di take on a finite number of values. (We call this set 
of possible values 3C). We show that the generating function of the sequence M(n,3,3*) is D-finite under 
certain conditions. 

A univariate generating function is D-finite if it satisfies a linear differential equation with polynomial 
coefficients. This property has considerable implications on asymptotic enumeration: D-finite functions 
have restrictions on their asymptotic form; asymptotic information is encoded in the differential equation 
that it satisfies; they can be treated with a number of automated tools. Although we do not compute 
asymptotic formulas here, the fact that the generating functions are D-finite can be useful for precisely such 
a computation. 

We cast this problem in graph theoretic language as follows in order to state our results precisely. Let 3 
and X be sets of positive integers, and suppose additionally that X is finite. Let Sg,x,n be the set of well- 
labelled graphs on n vertices where edge weights are from 3U{0} and the sum of weights of the edges incident 
to any given vertex is an element of X. The case when all di = k is the case of fc-regular graphs. Here a 
graph is well-labelled if the label set is { 1 , 2, ..., n}, where n is the number of vertices. 

The following is our main result. It appears below as Theorem [8] 

Theorem. Let 3 and X be finite sets of positive integers. Let Gg t %(z) be the generating function for the 
class 5g,x of well-labelled graphs where edge weights are from 3 U {0} and the sum of weights of the edges 
incident to any given vertex is an element of X. Then, 

Gg,x{z) = ^2 ISa,uc,rj| z n 

n 

is D-finite. Furthermore, the differential equation satisfied by Gg t x{z) is theoretically computable. 
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The condition of finiteness on 3 is not necessary, but is a consequence of the finiteness of X. We show 
how to relax the condition of finiteness, once there is more notation developed. 

For example G{i},{*,}(;?) is generating function for simple, labelled, fc-regular graphs, and G{ 1 2 ,...,fc},{fc}(- Z ) 
is the generating function for labelled fc-regular multigraphs. The D-finiteness of these generating functions 
for generic k was conjectured by Goulden, Jackson and Reilly [7], and was proved by Gessel [Bj. Our strategy 
proves the more general result by using the work of Mendez [13 to create the same framework as of [B], 
to which the work of [3] then applies. This amounts to building a symmetric function encoding of the full 
class of graphs, and then performing a subseries extraction to realize the degree restriction. These are all 
theoretically effective, and hence the differential equations are potentially computable. 

As Gessel noted, the class of all regular graphs is not D-finite. Consequently, our results, up to finite 
union of classes are most likely optimal. 

1.1. Contribution: Structure of graphs with controlled edge multiplicities. McKay remarked in his 
problem “The simplest non-trivial case is 3* = {0} and 3 = {0, 2, 3}”. In our notation, 3* = {0} corresponds 
to the criterion that the graphs have no loops, and our 3 are the same. We first consider di = k for all i, 
and loosen this to consider then di £ X for finite X and conclude with some more general comments. We 
also address 3* = {0,1}. 

Our principal contribution is a new formulation of his problem, and the resulting proof of D-finiteness. 
We find our approach to be of interest as it gives a new way to view McKay’s problem, and also a new 
example for the class of symmetric species. In this context, generalizations to hypergraph variants are very 
natural. 

The outline of the strategy has three main steps: 

(1) Find a series encoding a superset of graphs with edge weights from 3; 

(2) Rewrite this series as symmetric function using the power sum symmetric functions; 

(3) Extract the subseries of terms with the desired degree sequence via a scalar product operation. 

In Section [2] we describe how to write the graph generating functions using symmetric functions. The 
construction generalizes our previous work |14| in a straightforward way using the multiassemblies studied 
by Mendez m- The construction uses the species theory formalism [2], but we leave the category theoretic 
details to previous sources to avoid a rather substantial detour that is well described elsewhere. We remark, 
however, this appears to be different than the description of graphs as a species recently developed by Gainer- 
Dewar and Gessel [5j. That said, our Proposition [4] bears some resemblance to the formulas in Section 6 
of Henderson’s species generalization [8], which is the formalism upon which their work is based. Perhaps 
this problem is a good entry point for that theory. 

The D-finiteness result follows quickly once we adapt Waring’s formula. This is explained in Section [U 
We conclude with some directions on how to weaken the conditions as stated. 

2. Labelled graph generating functions 

We start with a systematic encoding of graph classes using symmetric functions. 

2.1. Simple graphs and X-generating functions. Let G be a simple graph with vertex set V{G) = 
{x\,... ,x n }, and edge set E(G). We associate to G the monomial i t(G) defined 

7r(G) = ]^[ XiXj = x^x^ 2 ■ ■. 

{®i,x^}€S(G) 

where di is the degree of x t . Let G(X) be the generating function of the set of all labelled simple graphs S, 
each with a vertex set a subset of X = {xi,X 2 

g(x)=5>(g)= n( i+ 

GeS i<j 

(Remark, this is a superset of the well-labelled graphs). To see the second equality, remark that every edge 
is either present once, or not at all. Similarly, if 9 is the set of graphs that permit multiple edges (but not 
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loops), 


g= v tt(g)= n 


1 


1 — XiX-j 


Ge 9 iK 'J 

as every edge exists some non-negative integer number of times. 

Under this description, the set S is a symmetric species [13], and the series encoding G(X) is what Mendez 
calls the associated X-generating function. Our strategy is to determine the X-generating function for the 
class of graphs of Sg in which edge weights are incorporated, and treated as multiplicities. (Hence, that the 
edge weights be positive integers is essential.) It is straightforward to get an expression for this, and then it 
is a mechanical manipulation to get a form we desire. 

The X-generating function for the class of labelled graphs with graphs edge weights from the set of positive 
integers 3 is 


( 1 ) 


Gg(X) = £ 7r(G) = 1 + E^) S 

GE Sg V sE8 


For example, G{ l! 2 , 3 ,...}(X) = G(X). When the X is clear, we might choose to not indicate the variable set. 

To determine the series for graph classes with loops it is sufficient to change the product index from i < j 
to i < j. If, rather, the loops have weights from a different set, say 3*, we multiply Eq. m by the product 
Eli (l + Sseg* ( x i x i) s ) ■ 

The number of well-labelled graphs in this class with given degree sequence d \,..., d n is simply the 
coefficient of the monomial x^x^ 2 ■ ■-x^ n . For suitably defined graph classes, and degree vectors, this is the 
value under investigation by McKay. In standard generating function notation we write the coefficient as 

[xtx d 2 *...x d n "]G g . 

As is always the case for symmetric species, the X-generating functions are symmetric functions. We access 
the coefficient using classic symmetric function operations. 


2.2. Expressing Gg(x) using symmetric functions. Since we can relabel any graph with a different set 
of labels and it remains in the class, Sg is a symmetric class with respect to the graph labels. We leverage 
this underlying symmetry to rewrite the generating function in terms of symmetric functions. 

Let A = (Ai,..., A*,) be a integer partition of n, a fact which we denote by A F n. Let X = xi, X 2 , ■ ■ ■ be 
an infinite, but countable, variable set. Then the symmetric function m\(X) (or simply m\) is defined 

m\{X) = E x ii x i2 ■ ■ ■ X ik 

where the sum is over all k- tuples of distinct positive integers (*i, * 2 , - - • 5 *fe)■ This is the monomial symmetric 
function indexed by A. The set of the monomial symmetric functions form a basis for a vector space of 
symmetric functions over Q. 

We express the classic elementary, complete, and power sum symmetric functions in the monomial basis 
as follows: 

h n — ^ [ 77 Za, Pn 771 (n) • 

Ah n 

Recall the definition e\ = e\ 1 e\ 2 ... e\ k . The set of the elementary symmetric functions indexed by partitions 
also forms a basis for a vector space of symmetric functions. This is also true for the sets of h\ and p\ 
respectively, which are similarly defined. We work in the ring of symmetric series over X: 

A((X)) := Q[[pi,p 2 ,P3, ■•■]]• 

In particular we are interested in symmetric series R(X) of the form 

r(a>EE c ^( x )- 

0<n Ahra 

The symmetric function operation known as plethysm is essential to our solution. Given two symmetric 
functions u and v, the inner law defines the quantity it [a] by defining the following rules, with u,v,w £ A 
and a, /3 in K 

(au + j3v)[w\ = au[w] + f3v[w], (-itu)[w] = Tt[w]u[w], 


3 



and, most importantly, if w = Y^,\°\P\ then p n [w\ = )T) A CAP( ra Ai)P(nA 2 ) • ■ For example, we can deduce 
that w[p n ] = p n [w ], and in particular that p n [pm\ = Pnm■ In a mnemonic way: 

(2) w\p n ] =w(p ln ,p2n,---,Pkn,---) whenever w = w(pi,p 2 , ... ,Pk, ■ ■ •)• 

Let H = ]T n h n and E = ]T) n e n - Gessel noted that G and G can both be expressed using pletliysm: 

(3) G = 1~[ (1+XjXj) = y^e„[e 2 ] = E[e 2 ]. 

i<j n 

(4) G = n(iz^:)=EMe»] = HN. 

i<j x 1 3/ n 

Given Eq. often plethysm expressions are easier to manipulate when the symmetric functions are 
written in the power sum basis. We do this next in Section [2.31 and this is followed by a discussion on how 
to derive the plethym expressions for Gg(A) in [3] 


2.3. Expressions in the power sum basis. We recall the following classic lemma as it guides our work. 
It shows how to express an infinite sum of h n as a function of power sum symmetric functions. 


Lemma 1 (Waring formula). The following equations are true: 

H = E h ” = II = ex P ) ’ 

n i<j 1 \ 0 <k / 

E = E en = IIo +x? -) = exp (E^ _1 ) fc+lpfe / fc ) • 

n i<j \0 <k / 

Proof. The proof is elementary series manipulations: 

log n = ^ log wVr = E E :r ^ k = E E x ^ k = E p^ k - 

0<4 1 0<4 1 0<4 0</c 0<fc 0<4 0 <k 


□ 


Indeed, the plethysms are easier to analyse given this form: 



= exp(^p fe [e 2 ]/fc) 

= eX P(E ^PklPl/ 2 -P2/2]) 

= eX P(E -P™))- 

We can similarly express E[e 2 ] (and indeed H[ft. 2 ], H[e 2 ]). 

3. Labelled graphs as a symmetric species 

One of the innovations of species theory I0E3, is a rigorous combinatorial interpretation of the plethysm 
operation in terms of natural compositions of combinatorial structures pQ. Plethysm as an analytic analog 
to composition has been well studied since Polya’s composition theorem. Asymmetric series of Labelle 03 
is also an important relative to the A'-generating functions that we seek. 

The combinatorial understanding of the composition in the particular case of H and E is formally devel¬ 
oped by Mendez [T3], and we gave a direct interpretation in the case of graphs, and other variants in fl4l . 
In particular, [14] contains is a description of graphs as multisort, and ultimately symmetric, species. The 
interpretation is as follows: a simple labelled graph is a set of edges. Edges are sets of atomic structures. 
Each atom is coloured a colour from the infinite set X = {a:i,a; 2 ...,}, and atoms of the same colour are 
identified to form a vertex. The combinatorial composition of a set of edges is reflected in A'-generating 
function by a plethysm. 
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In the notation of classic species, the class of all labelled multigraphs with loops allowed given by the 
multisort species E[E 2 ](X), and the associated cycle index series is H[/i 2 ]. Here E is the species of sets, 
and E 2 is the species of a set of cardinality two. Thus, a graph is a set of edges, but repetitions are 
not controlled in the multisort version. The multiassembly construction of Mendez controls repetition of 
elements. A multiassembly of type A, denoted M\, is a multiset where the multiplicities of the elements 
are prescribed by the parts of the partition A. For example, a multiassembly of type A = (1,1,..., 1) is a 
usual set without repetitions. Example 3.7 and Proposition 3.9 of m describe how to get the A'-generating 
function associated a composition of a multiassembly and a some object. In particular, the A-generating 
function of a multiassembly M A (X) is m A (X), and the composition is realized by plethysm. 

Example 2. We view graphs as multiassemblies of edges, where edge multiplicities are given by A. The edges 
themselves are multiassemblies of vertices. For example, the species M( 3 j2 )[M( 1i 1 )(A')] is the set of graphs 
with one edge of multiplicity two, and one edge of multiplicity three, under all possible positive integer 
labellings. There are two possible shapes, both under all possible valid labellings: The first graph is on three 
vertices, and is path. Each possible labelling contributes a monomial of the form x^x^xj, to the X-generating 
function. The second shape is a graph with four vertices, and two edges, and the graph is not connected. 
Each graph of this type contributes a monomial x\xjx\x\. The generating function is thus 7715 , 3,2 + 7713 , 3 , 2 , 2 - 
This is precisely 771 ( 3 , 2 ) [t71(i,i)]- 

The X-generating of Gg of S 3 requires a sum over all possible partitions with parts from 3- Towards a 
more compact notation, define the symmetric function /g, n as follows: 

(5) fs,n = Y TOA ‘ 

Ahn,;Ai£3 

To unravel the definition, remark /{ 2 , 3},6 = 77i 2 , 2 ,2 + 7773,3 and further note that /{i}, n = e n and /{i,..., n },„ = 
h n . We summarize these results in the following proposition. 

Proposition 3. Fix 3, a nonempty integer partition, and X = xi,x 2 ,..., and infinite, but countable set of 
labels. 

(1) The symmetric species of simple graphs S 3 is isomorphic to U A |- n;Ai e gM A [M (1 , 1) (X)]. 

(2) For any degree sequence d = (di, ..., d n ) with sum D = d\ + • ■ • + d n , the coefficient 

[xf-x^.-.x^] Y m A[e 2 ] 

AhD;Ai£3 

is precisely the number of graphs in Sg,N,n with degree sequence d. Equivalently, this is the coefficient 
of md when expanded in the multinomial basis of symmetric functions Q 

Similarly, the analogous class of graphs permitting loopH, denoted Sg, is isomorphic as species to 

u M a [M (1 , 1} (X) + M (2) (X)], 

Ahn;Aje3 

and the coefficient 

[xf'x^ ...X??] Y 777a[^7 2 ] 

XhD-XiG3 

is precisely the number of graphs in Sg N n with degree sequence d. 

3.1. Series expressions. To prove the D-finiteness results, we need the X-generating functions in a different 
format. We define Fg = /g, n and hence Fpi j = E and F{ 1j2 , } = H. Now, Fg — n 0<i (l + Ssea *<)• 
We can generalize Lemma [l] to express Fg in the power sum symmetric function basis. The proof follows 
from very basic manipulations. 


1 Willi d normalized to be decreasing if necessary. 

2 Graphs with loops of weight 1 corresponds to 3* = {0,1} in the original problem. 
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Proposition 4. Let S = {j i,... ,je} be a set of distinct positive integers. Then, 


F a 



where a n is the following sum taken over all compositions a = ( 01 , 0 : 2 ,...) of n such that each part is 
contained in S'- 


( 6 ) 


a n 


E 



parts(o) 


Here parts(o) is the number of parts in the composition. 


Proof. First we note that as the elements of S are all positive, a n is well defined since the number of such 
compositions is finite. Next we apply the same log-exp expansion, and some very basic coefficient extraction 
formulas: 


lQ g E 

0 <n 


bgjla+Eo 

0<2 s £3 


-E lo § 

0<i 


1 

1 -Esea(- a; f) 


EE 

0<i 0<k 


(_!)*+! 

k 


(E<) fc 

sea 


E E 

0 <i 0 <n 


a n p n . 

0 <n 


□ 


There are a few simplifications to note. If S — {1}, then a n = since there is only one term in the 
summation in Equation ([6|). When S = {1, 2,... }, a n = ^ as the sum is over all compositions, and we invoke 
a Mobius inversion argument. Further, when S = {si,..., s^} is finite, we can express this as follows: 


(7) 


0 <n 


fs,n = exp 



*1 12 


ll 


Ps 1 ii +S2*2H-b 


Theorem 5. Let S = {ji, ■ ■ ■ ,j(} be a finite set of £ distinct positive integers. Then Gg, the X-generating 
function for the symmetric species of labelled simple graphs with edge weights from S satisfies: 



with a n as defined in Proposition [3] 

Proof. This follows from Propositions [3] and 01 and the fact that e 2 = pf — P 2 - D 


Returning to our example: 
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4. D-finite symmetric series 


Recall that a series S £ A'[[xi,..., x n ]] is D-finite in xi,..., x n when the set of all partial derivatives and 
their iterates, <9 iH Hn F/dx l f ■ ■ ■ dx z ™ , spans a finite-dimensional vector space over the field AT(xi,..., x n ). 
This was generalized to an infinite number of variables by Gessel [6], who had symmetric functions in mind. 
A series S £ K[[xi,X 2 ,... ]] is D-finite in the x* if the specialization to 0 of all but a finite (arbritrary) choice 
of the variable set results in a D-finite function (in the finite sense). In this case, many of the properties of 
the finite multivariate case hold true. One notable exception is closure under algebraic substitution, which 
requires additional hypotheses. 

The definition is then tailored to symmetric series by considering the algebra of symmetric series as 
generated over Q by the set of power sum symmetric functions {pi,P 2 , ■ ■ ■ }• A symmetric series is called 
D-finite when it is D-finite as a function of the pf s. The applicability of this definition will be apparent in 
a moment. 

The two prototypcial examples, H and E are easily seen to be D-finite, as any such specialization of 
variables results in an exponential of a polynomial, which is clearly D-finite. Similarly, from the expression 
in Proposition |T] we see that the same argument will hold for any Fg. 

Theorem 6. For any set of positive integers 3, the series Fg and Gg are both D-finite symmetric series 
with respect to the p-basis. 

Proof. For both of these symmetric series any specialization of the p variables so that only a finite number 
are non-zero leaves an exponential of a polynomial, which is easily shown to be D-finite in the remaining 
variables. We immediately conclude the D-finiteness of both Fg and Gg, given the two previous results □ 


4.1. Extracting the generating functions. In the power notation for integer partitions, A = l ni .. .k nk 
indicates that i occurs n,; times in A, for i = 1,2,... ,k. The normalization constant 

2a := l in ni! ■ • • k nk Uk\ 

plays the role of the square of a norm of p\ in the following important formula: 

(8) (p\,P») = 5\^z x , 

where is 1 if A = p and 0 otherwise. 

The scalar product is useful for coefficient extraction because (toa,= <$a,/x- If we write S in the 
form Y^\ C\m\, then the coefficient of x^ 1 ... x£ fc in S is ca = ( S , h\). 

The closure under Hadamard product of D-finite series mi yields the consequence: 

Theorem 7 (Gessel [6]). Let S and T be elements ofQ[z][[pi,p 2 ,...]], D-finite in the pi’s and z, and further 
suppose that T involves only finitely many of the pi’s. Then (S, T) is D-finite as a function of z, provided 
it is well-defined as a power series. 

From this and Lemma [6] our main theorem follows almost immediately. 

Theorem 8. Let 3 and X be sets of positive integers, and suppose additionally that X is finite. Let Gg^%(z) 
be the generating function for the class Sg,ae of well-labelled graphs where edge weights are from 3 U {0} and 
the sum of weights of the edges incident to any given vertex is an element of X. Then, 

Gg,x(z) = ^2 |Sj,3C,n| z n 

n 

is D-finite. Furthermore, the differential equation satisfied by Gg t %{z) is theoretically computable. 


Proof. The series are combinatorial generating functions, and so they exist. We remark that 


(9) 




The first argument to the scalar product is D-finite by Lemma [6] Furthermore, hkZ k is a polynomial 

in the power sum basis and 2 when X is finite. Hence the second argument is rational, and D-finite. The 
result follows by Theorem [7] We address the computability in the next section. □ 
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Example 9. The series can be expanded: 

g { 3 , 2 } = F { 3 i2 }[e 2 ] = (m 2 + m 3 + m 2 , 2 + m 3i2 + 7712,2,2 + m 3i3 + ...) [e 2 ] 

= m 2j2 + m 3i3 + m.4,2,2 + 3 m 2 , 2 , 2,2 + m 5i3j2 + m 3 , 3 , 2 ,2 

+ f»6,2,2,2 + 2 7714,4,2,2 + 6 7714,2,2,2,2 + ^4,4,4 + ^6,3,3 + 15 7712,2,2,2,2,2 + 3 77l 3 , 3 , 3 , 3 + . . . 


To extract the generating function of 2-regular graphs examine the coefficients of 77i( 2 ,2,...,2) : 

G { 3 , 2 },{ 2 }(~) = / G { 3 , 2 }) = z 2 + 3z 4 + 15z 6 + 105z 8 + .... 

\ 0 <n / 

In this case, it is corresponds the number of matchings. Similarly, 

G {3,2},{2,3}(^0 = / G {3,2}) ^3^ 2 fc ~ / = 2z 2 + 7z 4 + 36z 6 + 429z 8 + .... 


0<n k— 0 


An extraction for any degree sequence is possible. The D-finiteness result can be generalized to handle 
infinite % provided Y^ n ^Cai-m-a ex' h\t n is a D-finite symmetric series. In general, to determine the generating 
function of graphs with a fixed set of degree sequences 7), it suffices to consider the series ^d~ parts ^ d) - 

We could also mark something other than number of vertices. As we noted, the symmetric series h n z n 
is D-finite, and in fact, for any finite fc, the series h^z n is D-finite. This would extract the generating 
function for the subclass of all regular graphs on k vertices from a given graph class with 17 marking the 
regularity. The resulting generating function is also D-finite. 


4.2. Comments on effective computation. There are two computational tools at hand to compute 
Gg,jc(z). One could iteratively expand /g,n[e 2 ] in the monomial basis as we did in the previous example. It 
might be slightly more efficient to expand the exponential expression. In practice, we were able to get some 
small results with this strategy. 

Alternatively, we can make use of the fact that Gessel’s result is effective: given the system of differential 
equations satisfied by symmetric series F and G, there are algorithms [3] to compute the differential equation 
satisfied by the scalar product. (Of course, at least one of F and G must contain other variables). It is 
straightforward to define the system satisfied by Gg , since it is expressed as an exponential of a polynomial. 
Consequently in theory we can compute the differential equation satisfied by G^,yc(z). In practice, using cur¬ 
rent algorithms, the computations are too resource intensive to deliver results when the number of variables 
is more than 5. We were able to confirm the correctness in small cases. 


5. Other generalizations 

By playing with the inner series in the plethysm, we can enumeration other families of objects, such as 
hyper graphs, or cyclic coverings of sets. The details are essentially given in [14; . 

As mentioned above, we could be interested in other kinds of functions for the growth of d. Other options 
should be possible to extract, provided the generating function for the extractor is D-finite. 

A different future direction would be to try to adapt the approach of de Panaheu and Ramos [¥j for 
multigraphs to the weighted edge versions. 
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